home *** CD-ROM | disk | FTP | other *** search
/ Magnum One / Magnum One (Mid-American Digital) (Disc Manufacturing).iso / d18 / nrpas13.arc / RZEXTR.DEM < prev    next >
Text File  |  1991-05-01  |  1KB  |  45 lines

  1. PROGRAM d15r7(input,output);
  2. (* driver for routine RZEXTR *)
  3. CONST
  4.    glimax=11;
  5.    glnmax=10;
  6.    glncol=7;
  7.    nv=4;
  8.    nuse=5;
  9. TYPE
  10.    glyarray = ARRAY [1..nv] OF real;
  11. VAR
  12.    glx : ARRAY [1..glimax] OF real;
  13.    gld : ARRAY [1..glnmax,1..glncol] OF real;
  14.    dum,xest : real;
  15.    i,iest,j : integer;
  16.    dy,yest,yz : glyarray;
  17.  
  18. (*$I MODFILE.PAS *)
  19. (*$I RZEXTR.PAS *)
  20.  
  21. BEGIN
  22. (* feed values from a rational function *)
  23. (* fn(x) := (1-x+x**3)/(x+1)**n *)
  24.    FOR i := 1 to 10 DO BEGIN
  25.       iest := i;
  26.       xest := 1.0/i;
  27.       dum := 1.0-xest+xest*xest*xest;
  28.       FOR j := 1 to nv DO BEGIN
  29.          dum := dum/(xest+1.0);
  30.          yest[j] := dum
  31.       END;
  32.       rzextr(iest,xest,yest,yz,dy,nv,nuse);
  33.       writeln;
  34.       writeln('iest  :=  ',i:2,'   xest =',xest:8:4);
  35.       write('Extrap. function: ');
  36.       FOR j := 1 to nv DO write(yz[j]:12:6);
  37.       writeln;
  38.       write('Estimated error:  ');
  39.       FOR j := 1 to nv DO write(dy[j]:12:6);
  40.       writeln
  41.    END;
  42.    writeln;
  43.    writeln('Actual values:    ',1.0:12:6,1.0:12:6,1.0:12:6,1.0:12:6)
  44. END.
  45.